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INTRODUCTION 

The  author  has  been  actively  engaged  in  applied  research,  design, 
manufacture,  and  testing  of  parachutes  and  parachute  systems  since 
f 1952.  In  this  period  many  improvements  have  been  accomplished  in  the 

prediction  of  parachute  opening  shock  performance  and  trajectory 
analysis.  The  coming  of  age  of  the  computer  was  a blessing  in  that 
it  permitted  detailed  modeling  of  parachute  deployments,  but,  at  the 
same  time,  in  many  cases,  it  reduced  the  need  to  understand  some  of 
the  underlying  basics  of  the  parachute  deployment  process.  One  of 
those  basics  is  a scale  factor  parameter.  The  development  of  the 
scale  factor  element  is  fundamental  to  the  understanding  of  parachute 
performance.  This  paper  is  a presentation  of  the  author's  current 
ideas  and  reasoning  on  what  is  required  to  permit  establishment  of  a 
means  of  scaling  parachute  deployments  and  trajectories.  A discussion 
of  what  variables  are  required  and  examples  of  existing  application 
and  proposed  method  of  application  are  offered.  This  will  be  accom- 
plished by  first  developing  the  mass  ratio  quantity  by  use  of  a 
horizontal  trajectory  with  a fully  deployed  parachute.  The  concept 
will  then  be  extended  to  show  applicability  to  inflating  parachutes, 
and,  finally,  methods  for  including  the  mass  ratio  as  a variable  in 
the  canopy  inflation  process  are  suggested.  This  paper  is  presented 
with  the  purpose  that  the  thoughts  expressed  may  assist  and  stimulate 
other  experimentors  in  the  development  of  future  advanced  analysis 
techniques . 

DEVELOPMENT  OF  THE  MASS  RATIO  CONCEPT 

In  medel  testing  of  airplanes,  missiles,  ships,  etc.  in  wind 
tunnels  and  towing  basins,  Reynold's  number,  Froude  number,  and  Mach 
number  are  accepted  scale  factors.  These  parameters,  together  with 
ballistic  coefficient  (W/CdSq)  and  surface  loading  (W/S0),  have  not 
provided  a general  correlation  of  parachute  test  data.  A viable 
scaling  parameter  for  general  use  in  deployable  decelerator  testing 
has  not  been  developed,  although  some  experimentors  have  successfully 
modeled  particular  cases.  A look  at  the  methods  of  testing  gives  us 
a clue  as  to  why  the  aforementioned  scale  factors  do  not  apply  to 
parachute  deployment.  In  the  testing  of  airplanes,  missiles,  ships, 
etc.,  a rigid  model  is  constructed  and  mounted  in  the  wind  tunnel, 
towing  basin,  etc.,  and  data  recorded  at  one  or  more  constant 
velocities.  Thus,  three  parameters  which  are  important  to  parachute 
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testing  are  eliminated.  First,  the  geometry  of  a deploying  parachute 
undergoes  dramatic  changes  as  compared  to  the  constant  geometry  of 
the  rigid  model.  Second,  the  deploying  parachute  geometry  is  a func- 
tion of  deployment  time.  Third,  the  velocity  profile  obtained  during 
a parachute  deployment  is  dependent  upon  the  mass  of  the  assembly 
being  retarded.  Transient  geometry,  time,  and  weight  are  not  factors 
in  the  aforementioned  scale  parameters.  This  leads  to  some  very 
important  conclusions: 

a.  Reynold's  number,  Froude  number,  and  Mach  number  are  valid 
scale  factors  for  airplanes,  missiles,  ships,  etc.,  because  they 
contain  the  variables  which  affect  performance. 

b.  These  established  scale  parameters  are  unsuitable  for  para- 
chute deployments  because  they  do  not  contain  important  variables 
which  affect  parachute  performance. 

c.  For  any  quantity  to  be  a valid  scale  parameter  in  any 
process,  it  must  contain  those  variables  which  affect  system  per- 
formance. 

An  example  of  a "mass  ratio"  ratio  scale  factor  can  be  developed 
in  a theoretical  horizontal  point  mass  trajectory.  For  simplicity, 
only  the  steady  state  (constant  drag  area,  CDSQ)  will  be  considered; 
however,  it  will  later  be  demonstrated  that  the  developed  mass 
ratio  is  also  applicable  for  parachute  deployments. 

Example  1:  Determine  the  velocity  profile  of  an  automobile 
which  is  being  retarded  on  a horizontal  road  by  a fully  deployed 
parachute.  Assume  negligible  automobile  aerodynamic  drag  and  road 
friction  forces. 

With  reference  to  Figure  1, 
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Integrating  and  solving  for  the  velocity  ratio 
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The  quantity  2W/pgVstCDSo  is  in  effect  a mass  ratio.  It  can  be 
visualized  as  in  Figure  2 to  be  the  ratio  of  the  mass  of  the  system, 
W/g,  to  the  mass  of  air  contained  in  a right  circular  cylinder  of 
face  area  CDS0,  length  V§  t,  and  mass  density  p.  If  this  quantity  is 
denoted  by  M,  then  Equation  (1)  becomes 


FIGURE  2 VISUALIZATION  OF  THE  MASS  RATIO  CONCEPT 


Figure  3 illustrates  the  concept  of  how  the  air  mass  is  affected  as 
the  automobile  moves  along  the  retarded  trajectory.  As  time 
increases,  additional  air  mass  is  affected,  and  at  each  instant,  a 
definite  velocity  ratio  and  mass  ratio  exist,  as  shown  in  Equation 
(2) . Note  that  the  length  Vst  in  Figure  3 is  not  a true  trajectory 
distance  because  of  the  use  of  the  constant  initial  velocity  Vs  for 
all  times.  This  mass  ratio  contains  the  basic  variables  (altitude, 
p;  parachute  aerodynamic  size,  CDS0;  system  mass,  W/g;  velocity,  Vs; 
and  time)  necessary  to  define  a scale  factor  for  deployable 
decelerator  application.  Apparent  mass  is  another  effect  which 
should  probably  be  included,  but  is  not  presently  sufficiently  under- 
stood or  defined  as  to  permit  inclusion  in  the  analysis. 
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As  the  automobile  is  retarded,  the  Reynold's,  Froude,  and  Mach 
numbers  vary  along  the  trajectory  as  the  velocity  ratio  changes. 

For  the  steady-state  trajectory,  only  the  ballistic  coefficient, 
W/C^Sq,  remains  constant  throughout  the  trajectory  and  can  be  used 
as  indicative  of  performance.  If  Reynold's,  Froude,  or  Mach  num- 
bers are  to  be  considered,  the  question  becomes  which  of  the  many 
numbers  should  be  used.  From  Equation  (2)  the  velocity  ratio  can 
be  determined  for  any  time,  t,  into  the  trajectory.  Once  the 
velocity  is  known,  the  Reynold's,  Froude,  and  Mach  numbers  can  be 
determined.  However,  before  these  three  scale  factors  could  be 
computed,  the  mass  ratio  had  to  be  calculated  first  in  order  to 
obtain  the  velocity  ratio.  Therefore,  the  mass  ratio,  even  though 
it  also  varies  along  the  trajectory,  once  again  is  the  superior 
scale  factor  because  it  directly  affects  performance  rather  than 
being  a secondary  calculation. 

Figure  4 illustrates  the  variation  of  the  various  scale  factors 
as  a function  of  time  for  the  automobile  of  example  (1)  with  the 
following  parameter  values: 

Vs  = 300  fps,  CDS0  = 25  ft2,  W = 500  lb,  p = 0.002378  slugs/ft3. 
The  Reynold's,  Froude,  and  Mach  numbers  are  defined: 


★ 

Reynold's  number 

Froude  number* 
Mach  number* 

The  Mass  Ratio 


pVL 


F-‘^r 


M_  - — 


M 


2W 


P«VStCDSo 


(3) 

(4) 

(5) 

(6) 


Some  of  the  mass  ratio  variables  appear  in  each  of  the  preceding 
scale  factors.  Then  it  is  possible  to  express  the  mass  ratio  M as 
some  function  of  (Re,Fr,Mn) . The  only  limitation  is  that  a single 
value  of  V=VS  be  used. 


From  Equation  (3) 


* The  characteristic  length,  L,  in  Reynold's  and  Froude  numbers  is 
taken  as  the  square  root  of  the  drag  area.  For  a discussion  of 
characteristic  length,  see  Appendix  A. 
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From  Equation  (4) 


Therefore , 


pv, 


s bmF2R. 


From  Equation  (5)  Vs  - Mne 


Substituting  these  quantities  in  the  mass  ratio,  Equation  (6! 


M 


2Wca 


M? 


g2t/i  CqS.  R.F2 


(7) 


Equation  (7)  does  not  appear  to  be  a convenient  way  to  express  the 
mass  ratio  scale  factor,  even  if  we  knew  the  instantaneous  values 
of  Re,  F r , and  Mn  as  a function  of  the  time,  t,  but  it  does  indi- 
cate that  modeling  parachute  performance  is  a complex  function  of 
Rg , Fr , and  Mj^  • 


MASS  RATIO  APPLICATION  TO  DEPLOYING  PARACHUTES 

The  term  "ballistic  coefficient"  is  defined  as  the  ratio  of  the 
system  weight  to  the  steady-state  drag  area  of  the  parachute 
(W/CqSq) . This  parameter  is  currently  used  as  a measure  of  the 
state  of  parachute  operation.  Ballistic  coefficient  values  greater 
than  approximately  40  are  considered  as  infinite  mass,  while  values 
of  approximately  0.3  are  finite  mass.  As  the  parachute  inflates, 
the  drag  area  varies  from  instant  to  instant.  If  the  instantaneous 
drag  area  can  be  determined  as  a function  of  deployment  time,  then 
the  variation  of  ballistic  coefficient  during  deployment  can  be 
defined. 

Appendix  B develops  an  empirical  deployment  drag  area  ratio 
expression  for  specified  solid  cloth  types  of  parachutes  as  a 
function  of  deployment  time  ratio,  which  is  expressed  in  Equation 
(8)  . 
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The  development  of  this  expression  is  a very  important  concept 
because  it  demonstrates  that  the  parachute  geometry  variation 
during  deployment  is  independent  of  altitude  and  velocity  and  is 
only  dependent  upon  the  initial  geometry  (n)  and  the  deployment 
time  ratio.  Berndt,  in  Reference  1,  confirmed  this  proposition  by 
testing  parachute  systems  at  various  altitudes  and  velocities. 


For  the  convenience  of  the  reader,  the  paper  in  which  Equation 
(8)  was  developed  is  included  as  Appendix  B.  A guide  for  the  use 
of  Appendix  B is  included  as  Appendix  C.  When  the  instantaneous 
drag  area  (CDS)  in  the  transient  ballistic  coefficient  (B.C.), 
Equation  (9),  is  replaced  by  Equation  (8),  the  variation  of  the 
ballistic  coefficient  during  deployment  is  shown  to  be: 


B. C.  - *- 
CDS 

(9) 

w 

/ ♦ "|2 

cds„ 

(10) 

When  the  initial  condition,  n>  approaches  zero,  the  ballistic 
coefficient  approaches  infinity  at  time  t=o.  From  this  initial 
value,  the  ballistic  coefficient  then  modifies  from  instant  to 
instant  to  whatever  final  ballistic  coefficient  (W/C^Sq) (infinite 
mass,  intermediate  mass,  or  finite  mass)  the  particular  system 
under  study  possesses.  For  any  given  value  of  W/CqSq  and  n*  the 
variation  during  deployment  is  always  the  same  (see  Fig.  5) . 

Appendix  B developed  a mass  ratio  scale  factor  similar  to 
Equation  (6)  for  an  inflating  solid  cloth  parachute,  which  was  a 
constant  value  throughout  the  deployment  process.  The  only 
difference  between  the  mass  ratio  of  Appendix  B and  the  mass  ratio 
visualized  in  Figure  2 is  that  the  variable  time  t is  replaced 
by  the  discreet  time  tQ. 

The  magnitude  of  the  mass  ratio  (M)  together  with  the  initial 
area  effect  (o)  determine  the  state  of  parachute  operation 


1.  Berndt,  R.J.  and  DeWeese,  J.H.,  "Filling  Time  Prediction 

Approach  for  Solid  Cloth  Type  Parachute  Canopies,"  A1AA  Aero- 
dynamic Decelerations  Systems  Conference,  Houston,  Texas, 

7-9  Sep  1966. 
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(infinite  mass,  etc.),  the  velocity  ratios  during  deployment,  the 
maximum  opening  shock  force,  and  the  time  of  occurrence  of  the 
maximum  shock  during  or  after  the  unfolding  phase  of  deployment. 
Examination  of  Figures  9 and  10,  p.  B-5,  Appendix  B,  discloses 
that  in  this  system  of  measurement  mass  ratios  greater  than  10  are 
essentially  infinite  mass  state  of  deployment,  and  for  n=0  mass 
ratios  less  than  0.19  are  finite  mass  state  of  deployment.  In 
between  these  limits  are  intermediate  mass  ratio  values  where  the 
parachute  opening  shock  occurs  after  time  tQ,  but  can  be 
significantly  less  than  infinite  mass  calculations.  The  mass 
ratio  limit  for  the  finite  state  of  deployment  is  dependent  on 
the  value  of  n*  An  analysis  of  limiting  mass  ratios  is  presented 
in  Appendix  D. 

In  the  method  of  developing  the  inflating  parachute  analysis, 
only  one  value  of  mass  ratio  was  possible.  The  concept  of  a non- 
constant mass  ratio  during  inflation  was  an  innovation  which 
occurred  after  the  publication  of  Appendix  B.  A possible  develop- 
ment of  variable  mass  ratio  during  canopy  inflation  requires  that 
the  reference  time  tQ  and  the  steady-state  drag  area  C^Sq  be 
replaced  by  their  variable  counterparts  t and  CUS. 


pgvstcDs 


When  CDS  is  defined  by  Equation  (8)  and  the  denominator  is  multi- 
plied By  t0/tQ=l. 


/t\3 

~]  2 

-V.V.g) 

+ T? 

(11) 


The  mass  ratio,  unlike  the  ballistic  coefficient,  is  infinite  for 
all  values  of  n at  the  beginning  of  deployment.  The  variations 
of  mass  ratio  during  deployment  is  compared  to  the  variation  of 
ballistic  coefficient  as  a function  of  n in  Figure  5. 

A comparison  of  the  development  of  the  affected  air  volumes 
for  single  value  and  variable  mass  ratios  is  visualized  in  Figure 
6.  The  variable  mass  ratio  air  volume,  as  shown,  is  for  a mass 
ratio  greater  than  the  limiting  value.  In  the  method  of  calcula- 
tion in  Appendix  B,  the  elasticity  of  the  materials  are  not  con- 
sidered significant  in  finite  mass  deployments,  because  the 
maximum  opening  shock  occurs  before  full  inflation. 
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CONCLUSIONS 


The  following  conclusions  are  derived  from  the  arguments  pre- 
sented in  this  paper. 

1.  Scale  factors  for  use  in  parachute  performance  analysis 
must  contain  those  variables  which  affect  operation.  Weight,  time, 
variable  geometry,  and  air  flow  through  the  canopy  are  very  impor- 
tant parameters  not  accounted  for  by  Reynold's,  Froude,  and  Mach 
number  similitude  expressions.  The  mass  ratio,  M,  which  does 
contain  the  basic  parameters  required  for  a meaningful  scale 
factor,  was  derived  from  the  analysis  of  a horizontal  retarded 
trajectory.  The  visualization  of  Equation  (6),  in  Figure  2,  when 
viewed  in  this  manner,  demonstrates  that  this  quantity  is  a true 
ratio  of  masses,  while  Equation  (2)  and  Figure  3 show  that  the 
trajectory  performance  is  dependent  upon  the  affected  air  mass  as 

a function  of  time. 

2.  Scale  factors  for  parachutes  are  likely  to  be  complicated 
functions  of  Reynold's,  Froude,  and  Mach  numbers.  An  improvement 
in  the  application  of  Reynold's  number  and  Froude  number  to  para- 
chutes can  be  attained  by  utilizing  a characteristic  length  whose 
magnitude  is  equal  to  the  square  root  of  the  steady-state  drag 
area.  This  reference  length  provides  a logical  basis  for  comparing 
the  retarded  fully  inflated  parachute  trajectories  for  any  type  of 
parachute  for  a single  value  of  Reynold's  number.  At  the  same 
time,  canopy  deployments  for  parachutes  of  differing  inflation 
characteristics  will  vary  for  identical  Reynold's  number  or  Froude 
number . 


New  names  were  applied  to  the  steady-state  drag  area  and 
the  characteristic  length.  These  names  are  "Aerodynamic  Size"  and 
"Aerodynamic  Diameter,"  respectively. 

3.  The  realization  that  the  geometry  of  parachute  deployment 
is  independent  of  altitude,  velocity,  and  system  mass  is  a very 
important  discovery!  The  deployment  geometry  can  now  be  separated 
from  the  forces  generated  and  analyzed  separately.  Drag  area 
signatures  obtained  from  infinite  mass  wind  tunnel  tests  are 
applicable  to  finite  mass  and  intermediate  mass  analysis. 

4.  The  mass  ratio  equation,  M,  is  the  same  for  all  types  of 
parachutes.  When  different  types  of  parachutes  are  tested  under 
identical  deployment  conditions,  the  particular  mass  ratio  for  a 
given  parachute  type  is  dependent  upon  the  reference  time  tQ. 

The  reference  time  tQ  is  dependent  upon  the  geometry,  air  flow 
properties  of  the  canopy,  drag  area  signature,  and  deployment 
conditions.  For  the  solid  cloth  types  of  parachutes  of  Appendix 
B,  the  dependency  of  the  mass  ratio  and  reference  time  on  the 
aforementioned  parameters  are  illustrated  in  the  flow  chart  of 
Figure  7.  The  technique  of  Appendix  B can  be  used  to  develop  an 
analysis  for  other  types  of  parachutes,  once  the  inflation  drag 
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FIGURE  7 DEPENDENCE  OF  MASS  RATIO  AND  REFERENCE  TIME 
ON  PARACHUTE  GEOMETRY,  AIR  FLOW  PROPERTIES, 
DRAG  AREA  SIGNATURE,  AND  DEPLOYMENT 
CONDITIONS  FOR  SOLID  CLOTH  PARACHUTES 
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area  signature  has  been  determined  as  a function  of  initial  area 
and  deployment  time  ratio  from  infinite  mass  wind  tunnel  deploy- 
ments and  methods  of  handling  airflow  through  the  canopy  are 
developed.  Treating  the  drag  area  as  a single  quantity  as  opposed 
to  instantaneous  drag  coefficients  and  areas  simplifies  the  prob- 
lem and  appears  to  be  adequate  for  calculating  forces. 

5.  The  mass  ratio  concept  from  the  horizontal  trajectory 
analysis  was  applied  to  the  inflation  analysis  of  several  types 
of  solid  cloth  parachutes  of  Appendix  B.  A single  value  of  mass 
ratio  was  used  throughout  the  inflation  process,  and  reasonable 
results  were  obtained.  In  addition,  the  effects  of  altitude  on 
inflation  distance  and  shock  factor  were  analytically  described. 

6.  As  a further  advance  in  inflation  analysis,  the  idea  of  a 
variable  mass  ratio  during  canopy  inflation  was  introduced.  For 
both  ballistic  coefficient  and  mass  ratio  approaches,  it  was 
demonstrated  that  all  inflation  processes  begin  as  infinite  mass 
quantities.  Future  analysis  of  opening  shock  which  utilizes  a 
variable  mass  ratio  during  deployment  should  improve  the  accuracy 
of  the  analysis  as  it  more  nearly  approaches  the  realistic 
environment . 

7.  The  initial  area  of  the  parachute  at  suspension  line 
stretch  is  dependent  upon  the  deployment  system  for  magnitude  and 
repeatability.  The  effects  of  n on  deployment  forces  can  be 
evaluated  by  treating  n as  a variable  and  calculating  the  result- 
ing forces  for  constant  deployment  conditions. 

As  n increases,  the  reference  time  tQ  is  reduced  with  a 
subsequent  rise  in  the  mass  ratio  and  shock  factor.  In  finite 
mass  deployments,  the  increase  in  parachute  drag  force  early  in 
inflation  reduces  the  maximum  force  at  t = tQ.  The  amount  of 
reduction  depends  upon  the  mass  ratio.  At  a value  of  approximately 
n = 0.4,  the  shock  factors  at  t = 0 and  t = tQ  are  essentially 
equal  (see  Fig.  D-2) . Initial  area  ratios  of  this  magnitude  can 
be  adapted  to  the  deployment  analysis  of  the  second  stage  of  a 
reefed  parachute  by  considering  the  disreefing  as  the  deployment 
of  a parachute  with  a large  n* 

8.  As  the  mass  ratio  increases  during  finite  deployment,  the 
maximum  force  occurs  later  in  the  inflation  process.  The  mass 
ratio  which  causes  the  maximum  force  to  occur  at  t = tQ  is 
defined  as  the  limiting  mass  ratio,  ML,  which  is  a function  of  n . 


21 


NSWC/WOL  TR  78-189 


Appendix  A 

DISCUSSION  OF  CHARACTERISTIC  LENGTH 


The  application  of  Reynold's  number  and  Froude  number  to  the 
scaling  of  parachute  performance  is  not  satisfactory.  One  cause 
of  this  deficiency  is  the  definition  of  the  characteristic  length 
used  in  each  of  these  parameters.  Many  experimenters  quote  values 
of  Reynold's  number  per  foot  of  length,  regardless  of  the 
geometric  size  of  the  parachute,  in  order  to  avoid  the  necessity 
of  defining  a characteristic  length. 

If  a wide  range  of  any  given  type  of  parachute  sizes  were 
tested  at  constant  deployment  conditions  and  attached  mass,  the 
opening  shock  characteristics  and  the  inflight  trajectories  vary 
for  each  of  the  parachute  sizes  used,  but  the  Reynold's  number  per 
foot  remains  constant.  This  is  illogical.  It  could  be  argued  in 
this  case  that  use  of  the  parachute  design  diameter,  D0,  would 
suffice  as  a characteristic  length  because  the  various  trajectories 
are  different  by  reason  of  the  different  canopy  sizes.  Where  this 
argument  breaks  down  is  when  a different  type  of  parachute  of  the 
same  design  drag  area  (CqSq)  is  substituted,  such  as  a ribbon  type 
of  canopy  being  substituted  for  a solid  flat  canopy.  In  this  case 
the  inflight  trajectories  are  matched,  but  the  geometric  diameter 
of  the  ribbon  chute  is  considerably  larger  than  the  equivalent 
flat  circular.  The  use  of  the  geometric  diameter  in  this  case 
would  give  two  values  of  Reynold's  number  for  the  same  inflight 
trajectory.  The  "best  choice"  of  reference  length  is  that  length 
which  will  define  the  inflight  trajectories  for  similar  ballistic 
coefficients  (W/CdSq)  by  one  value  of  Reynold’s  number.  The 
characteristic  length  described  in  Equation  (A-l)  meets  the  single 
value  criteria. 


L - D 


A 


(A-1) 


The  opening  shock  force  characteristics  of  the  different  types  of 
parachutes  will  vary  for  identical  Reynold's  numbers.  The  fact 
that  the  deployment  performance  of  various  types  of  parachutes  is 
different  for  identical  Reynold's  numbers  is  to  be  expected  and 
logical . 

The  following  proposals  are  offered: 


A-1 


NSWC/WOL  TR  78-189 


a.  The  characteristic  length  used* in  Reynold's  number  calcula- 
tions as  applied  to  parachutes  be  defined  as  per  Equation  (A-l)  and 
be  named  the  "Aerodynamic  Diameter."  The  aerodynamic  diameter  shall 
be  designated  by  the  symbol  "DA  » 

b.  If  we  consider  the  surface  area  (SQ)  of  a parachute  as 
related  to  the  geometric  size  of  the  canopy,  then  in  che  same 
manner  the  drag  area  (CqSq)  can  be  related  to  the  designation 
"Aerodynamic  Size." 


NSWC/WOL  TR  78-189 


► 


, 


4 


Appendix  B 


TECHNIQUE  FOR  THE  CALCULATION  OF  THE 
OPENING-SHOCK  FORCES  FOR  SEVERAL 
TYPES  OF  SOLID  CLOTH  PARACHUTES 


AIAA  Paper 
No.  73-477 


A TECHNIQUE  FOR  THE  CALCULATION  OF  THE 
OPENING-SHOCK  FORCES  FOR  SEVERAL  TYPES  OF 
SOLID  CLOTH  PARACHUTES 

by 

W.  P.  LUDTKE 

Naval  Ordnance  Laboratory 

Silver  Spring,  Maryland 


AIM  4m  Aerodynamic 
Deceleration  Systems 
Conference 

PALM  SPRINGS,  CALIFORNIA  / MAY  21-23,  7973 

First  publication  rights  reserved  by  American  Institute  of  Aeronautics  and  Astronautics. 

1290  Avenue  of  the  Americas.  New  York,  N.  Y.  10019.  Abstracts  may  be  published  without 
permission  if  credit  is  given  to  author  and  to  AIAA  (Price  AIAA  Member  SI. 50.  Nonmember  $2.00). 

Note:  This  paper  available  at  AIAA  New  York  office  for  si*  months, 
thereafter,  photoprint  copies  are  available  at  photocopy  prices  from 
AIAA  Library,  750  3rd  Avenue,  New  York,  New  York  10017 


NSWC/WOL  TR  78-189 


A TECHNIQUE  FOR  THE  CALCULATION  OF  THE  OPENING-SHOCK  FORCES  FOR 
SEVER.i  TYPES  OF  SOLID  CLOTH  PARACHUTES 

W.  P.  Ludtke 

Naval  Ordr.arce  Laboratory 
Silver  Spring,  Maryland 


r 


I 

I 


I 


Abstract 


An  analytical  method  of  calculating 
parachute  opening-shock  forces  based  upon 
wind-tunnel  derived  drag  area  time  signa- 
tures of  several  solid  cloth  parachute 
types  in  conjunction  with  a scale  factor 
and  retardation  system  steady-state  param- 
eters has  been  developed.  Methods  of 
analyzing  the  Inflation  time,  geometry, 
cloth  airflow  properties  and  materials 
elasticity  are  included.  The  effects  of 
mass  ratio  and  altitude  on  the  magnitude 
and  time  of  occurrence  of  the  maximum 
opening  shock  are  consistent  with  observed 
field  test  phenomena. 

I.  Introduction 


In  1965,  the  Naval  Ordnance  Laboratory 
(NOL)  was  engaged  in  a project  which 
utilized  a 35-foot-diameter,  10-percent 
extended-skirt  parachute  (type  T-10)  as 
the  second  stage  of  a retaraation  system 
for  a 250-pound  payload.  Deployment  of 
the  T-10  parachute  was  to  be  accomplished 
at  an  altitude  of  100,000  feet.  In  this 
rarefied  atmosphere,  the  problem  was  to 
determine  the  second  stage  deployment 
conditions  for  successful  operation.  A 
search  of  available  field  test  information 
indicated  a lack  of  data  on  the  use  of 
solid  cloth  parachutes  at  altitudes  above 
30,000  feet. 

The  approach  to  this  problem  was  as 
follows:  Utilizing  existing  wind-tunnel 
data,  low-altitude  field  test  data,  and 
reasonable  assumptions,  a unique  engi- 
neering approach  to  the  inflation  time 
and  opening-shock  problem  was  evolved  that 
provided  satisfactory  results.  Basically, 
the  method  combines  a wind-tunnel  derived 
drag  area  ratio  signature  as  a function 
of  deployment  time  with  a scale  factor  and 
Newton's  second  law  of  motion  to  analyze 
the  velocity  and  force  profiles  during 
deployment.  The  parachute  deployment 
sequence  is  divided  into  two  phases.  The 
first  phase,  called  "unfolding  phase," 
where  the  canopy  is  undergoing  changes  in 
shape,  is  considered  to  be  inelastic  as 
the  parachute  inflates  initially  to  its 
steady-state  aerodynamic  size  for  the 
first  time.  At  this  point,  the  "elastic 
phase"  is  entered  where  it  is  considered 
that  the  elasticity  of  the  parachute 
materials  enters  the  problem  ar.d  resists 
the  applied  forces  until  the  canopy  has 
reached  full  inflation. 

The  developed  equations  are  in  agree- 
ment with  the  observed  performance  of 
solid  cloth  parachutes  in  the  field,  such 
as  the  decrease  of  inflation  time  as 


altitude  increases,  effects  of  altitude  on 
opening-shock  force,  finite  and  Infinite 
mass  operation,  and  inflation  distance. 

II.  Development  of  Velocity  Ratio  and 
Force  Ratio  Equations  During 
the  Unfolding  Phase  of 
Parachute  Deployment 

The  parachute  deployment  would  take 
place  in.  a horizontal  attitude  in  accord- 
ance with  Newton's  second  law  of  motion. 

EF  = ma 


\ P y2cDS 


W dV 
g dt 


It  was  recognized  that  other  factors,  such 
as  included  air  mass,  apparent  mass,  and 
their  derivatives,  also  contribute  forces 
act!  ■»  on  the  system.  Since  definition 
of  t..ese  parameters  was  difficult,  the 
analysis  was  conducted  in  the  simplified 
form  shown  above.  Comparison  of  calculated 
results  and  test  results  indicated  that 
the  omitted  terms  have  a small  effect. 
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Multiplying  the  right-hand  side  of  equa- 
tion (1)  by 
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and  rearranging 
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In  order  to  integrate  the  left-hand  term 
of  equation  (2),  the  drag  area  ratio  must 
be  defined  for  the  type  of  parachute  under 
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analysis  as  a function  of  deployment 
reference  time,  tQ. 


TIME  RATIO  1 1<„ 

FIG.  1 TYPICAL  INFINITE  MASS  FORCE-TIME  HISTORY  OF  A 
SOLID  CLOTH  PARACHUTE  IN  A WIND  TUNNEL 


TIME  (SECI 

REPRODUCED  FROM  REFERENCE  ID 

FIG.  3 TYPICAL  FORCE  TIME  CURVE  FOR  A 10%  EXTENDED 

SKIRT  PARACHUTE  UNDER  INFINITE  MASS  CONDITIONS. 


FIG.  4 TYPICAL  FORCE  TIME  CURVE  FOR  A PERSONNEL 
GUIDE  SURFACE  PARACHUTE  UNDER  INFINITE 
MASS  CONDITIONS 


Figure  1 illustrates  a typical  solid 
cloth  parachute  wir.d-tunnel  infinite  mass 
force-time  history  after  snatch.  In 
infinite  mass  deployment,  the  maximum  size 
ar.d  maximum  shock  force  occur  at  the  time 
of  full  inflation,  tf.  However,  tf  is 
inappropriate  for  analysis  since  it  is 
dependent  upon  the  applied  load,  structural 
strength,  and  materials  elasticity.  The 
reference  time,  t,-,  where  the  parachute  has 
attained  its  steady-state  aerodynamic  size 
for  the  first  time,  is  used  as  the  basis  for 
performance  calculations. 

At  any  instant  during  the  unfolding 
phase,  the  force  ratio  F/Fs  can  be  deter- 
mined as  a function  of  the  time  ratio, 
tAc  • 


tpVs  CDSo 


Since  the  wind-tunnel  velocity  and  density 
are  constant  during  infinite  mass  deployment 


CDS 

CDSo 


TIME  (SECI 


FIG.  2 TYPICAL  FORCE  TIME  CURVE  FOR  A SOLID  FLAT  PARA 
CHUTE  UNDER  INFINITE  MASS  CONDITIONS 


FIG  5 TYPICAL  FORCE  TIME  SIGNATURE  FOR  THE  ELLIPTICAL 
PARACHUTE  UNDER  INFINITE  MASS  CONDITIONS 


FIG  6 TYPICAL  FORCE  TIME  SIGNATURE  FOR  THE  RING 

SLOT  PARACHUTE  20%  GEOMETRIC  POROSITY  UNDER 
INFINITE  MASS  CONDITIONS 


Infinite  mass  opening-shock  signatures 
c ' several  types  of  parachutes  are  pre- 
se-teJ  in  figures  2 through  6.  Analysis 
of  these  signatures  using  the  force  ratio, 
F/F-,  - tire  ratio,  tA0>  technique 
'-diea'ed  a similarity  in  the  performance 
of  the  various  solid  cloth  types  of 
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parachutes  which  were  examined.  The  geo- 
metrically porous  ring  slot  parachute  i 

displayed  a completely  different  signature, 
as  was  expected.  These  data  are  illus- 
trated in  Figure  7.  If  &n  initial  boundary 


EQ.  4,  ij  - 0 
- ea  «.  i)  - o.2 

PERSONNEL  GUIDE  SURFAC^j 
RING  SLOT  CANOPY  f 

. ELLIPTICAL  CANOPY  — - I 
> SOLID  FLAT  CANOPY  A 

10%  EXTENDED  SKIRT  / 
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FIG.  7 DRAG  AREA  RATIO  VS.  TIME  RATIO 

condition  of  C^S/C^S0  = 0 at  tine  t/tQ  = 0 

is  assumed,  then,  the  data  can  be  approxi- 
mated by  fitting  a curve  of  tne  form 


--  ‘ (“)6 

' DSo  Vo/ 
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At  the  tine  that  equation  (5)  was  ascer- 
tained, It  suggested  that  the  geometry  of 
the  deploying  parachute  was  Independent 
of  density  and  velocity.  It  was  also 
postulated  that  although  this  expression 
had  been  determined  for  the  infinite  mass 
condition,  it  would  also  be  true  for 
the  finite  mass  case.  Tnis  phenomenon 
has  since  been  independently  observed  and 
confirmed  by  Berndt  and  De  Weese  In 
reference  (2). 

Since  the  drag  area  ratio  was  determined 
from  actual  parachute  deployments,  it  was 
assumed  that  the  effects  of  apparent  mass 
and  included  mass  on  the  deployment  force 
history  were  accommodated. 

The  right-hand  term  of  equation  (2) 
contains  the  expression 


P8VoCDSo 


This  term  can  be  visualized  as  shown  in 
Figure  8 to  be  a ratio  of  the  retarded 
mass  (including  the  parachute)  to  an 
associated  mass  of  atmosphere  contained 
in  a right  circular  cylinder  which  is 
generated  by  moving  an  inflated  parachute 
of  area  CjjS0  for  a distance  equal  to  the 
product  of  Vst0  through  ar.  atmosphere  of 
density,  o. 


CYLINDER 
FACE  AREA 


VOLUME  DENSITY  ■ 


A mo*e  realistic  drag  area  ratio  expres- 
sion wan  determined  which  Includes  the 
effect  of  initial  area  at  line  stretch. 


where  n is  the  ratio  of  the  projected 
mouth  area  at  line  stretch  to  the  steady- 
state  projected  frontal  area.  Expanding 
equation  (4) 


2V» 
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FIG  8 VISUALIZATION  OF  THE  MASS  RATIO  CONCERT 


The  mass  ratio,  M,  is  the  scale  factor 
which  controls  the  velocity  and  force 
profiles  during  parachute  deployment. 
Substituting  M and  CdS/CdSq  into  equa- 
tion (2),  Integrating,  and  solving  for 
V/V„ 
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The  Instantaneous  shock  factor  Is  defined 

as 


x 


1 


fpygy 

i^sVo 


Figures  9 and  10  illustrate  the  velocity 
and  force  profiles  .generated  from  equations 
(7)  and  (8)  for  initial  projected  area 
ratios  of  n = 0,  and  0.2  with  various  mass 
ratios. 

IV.  Methods  for  Calculation  of 
the  Reference  Time,  tG 


If  the  altitude  variation  during  deployment 
is  small,  then,  the  density  may  be  consid- 
ered as  constant 


from  equations  (5)  and  (7) 


III.  Maximum  Shock  Force  and  Time  of 
Occurrence  During  the  Unfolding  Phase 

The  time  of  occurrence  of  the  maximum 
Instantaneous  shock  factor,  x^,  is  diffi- 
cult to  determine  for  the  general  case. 
However,  for  n = 0,  the  maximum  shock 
factor  and  time  of  occurrence  are  readily 
calculated.  For  n = 0 


The  ratio  concept  is  an  ideal  method  to 
analyze  the  effects  of  the  various  param- 
eters on  the  velocity  and  force  profiles 
of  the  opening  parachutes;  however,  a 
means  of  calculating  t0  is  required  before 
specific  values  can  be  computed.  Methods 
for  computing  the  varying  mass  flow  into 
the  inflating  canopy  mouth,  the  varying 
mass  flow  out  through  the  varying  inflated 
canopy  surface  area,  and  the  volume  of 
air,  V_,  which  must  be  collected  during 
the  Inflation  process  are  required. 

Figure  11  represents  a solid  cloth- 
type  parachute  canopy  at  some  instant 
during  Inflation.  At  any  given  instant, 
the  parachute  drag  area  is  proportional 
to  the  maximum  inflated  diameter.  Also, 
the  maximum  diameter  in  conjunction  with 
the  suspension  lines  determines  the  inflow 
mouth  area  (A-A)  and  the  pressurized 
canopy  area  (B-B-B).  This  observation 
provided  the  basis  for  the  following 
assumptions.  The  actual  canopy  shape  is 
of  minor  importance. 

a.  The  ratio  of  the  instantaneous 
mouth  inlet  area  to  the  steady-state 
mouth  area  is  in  the  same  ratio  as  the 
instantaneous  drag  area. 


Setting  the  derivative  of  x^  with  respect 
to  time  equal  to  zero  and  solving  for 
t/tQ  at  xt  max 


max 


(9) 


and  the  maximum  shock  factor  is 
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(10) 
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b.  The  ratio  of  the  Instantaneous 
pressurized  cloth  surface  area  to  the 
canopy  surface  area  is  in  the  same  ratio 
as  the  Instantaneous  drag  area. 


CpS_ 

CDSo 


c.  Since  the  suspension  lines  in  the 
unpressurized  area  of  the  canopy  are 
straight,  a pressure  differential  has  not 
developed,  and,  therefore,  the  net  air- 
flow in  this  zone  is  zero. 

Based  on  the  foregoing  assumptions, 
the  mass  flow  equation  can  be  viritten 


Equations  (9)  and  (10)  are  valid  for  valuta 

°f  M * iT  (0-19)»  8lnce  for  larger  values 

of  M,  the  maximum  shock  force  occurs  in 
the  elastic  phase  of  Inflation. 


dm  = m inflow  - m outflow 


P 77  = Pvam  - pAsp 


V/V,  VELOCITY  RATIO 
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\ 
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FIG.  11  PARTIALLY  INFLATED  PARACHUTE  CANOPY 


From  equation  (3) 

i°l . (±\6 
Vo  v*J 

From  equation  (7) 

V. 

V - 


; for  r|  = 0 


-(--Y 

™V‘o  / 


; n = o 


l + 


From  equation  (26) 


r2n 


Integrating: 


Measured  values  of  n Indicate  a data 
range  from  0.574  through  0.771 . A 
convenient  solution  to  the  reference  time 
equation  evolves  when  n Is  assigned  a 
value  of  1/2.  Integrating  equation  (13) 
and  using 


Equation  (14)  expresses  the  unfolding 
reference  time,  tQ,  In  terms  of  mass, 
altitude,  snatch  velocity,  airflow  char- 
acteristics of  the  cloth,  and  the  steady- 
state  parachute  geometry.  Note  that  the 
term  &pV_(/ W is  the  ratio  of  the  Included 
air  mass  to  the  mass  of  the  retarded 
hardware.  Multiplying  both  sides  of 
equation  (14)  by  Vs  demonstrates  that 


V t = a constant  which  is 

a function  of  altitude 


Figures  1?  and  13  Indicate  the  para- 
chute unfolding  time  and  unfolding 
distance  for  values  of  n = 1/2  and  n = 
0.63246.  Note  the  variation  and  conver- 
gence with  rising  altitude.  The  opening- 
shock  force  is  strongly  Influenced  by 
the  Inflation  time.  Because  of  this,  the 


i 
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value  of  t0  calculated  by  using  a real- 
istic value  of  n should  be  used  in  the 
lower  atmosphere. 

As  an  example  of  this  method  of  opening- 
shock  analysis,  let  us  examine  the  effect 
of  altitude  on  the  opening-shock  force  of 
a T-10-type  parachute  retarding  a 200- 
pound  weight  from  a snatch  velocity  of 
Vs  = 400  feet  per  second  at  sea  level. 
Conditions  of  constant  velocity  and  con- 
stant dynamic  pressure  are  investigated. 

The  results  are  presented  in  Figure  14. 

At  low  altitudes,  the  opening- shock  force 
is  less  than  the  steady-state  drag  force: 
however,  as  altitude  rises,  the  opening 
shock  eventually  exceeds  the  steady-state 
drag  force  at  some  altitude.  This  trend 
is  in  agreement  with  field  test  observations. 

V.  Correction  of  tD  for 
Initial  Area  Effects 


The  unfolding  reference  time,  tQ, 
calculated  by  the  previous  methods  assumes 
that  the  parachute  inflates  from  zero  drag 
area.  In  reality,  a parachute  has  a drag 


UNFOLDING  TIME  (SEC  I 


FIG.  12  EFFECT  OF  ALTITUDE  ON  THE  UNFOLOING  TIME 
~l  '•  AT  CONSTANT  VELOCITY  AND  CONSTANT 
DYNAMIC  PRESSURE  FOR  n - 1/2  AND  n * 0 63296 


CLOTH  MIL  C 7020.  TYPE  III 

3b  FOOT.  10%  EXTENOEO  SKIRT  CANOPY 


V.'o 


FIG  13  EFFECT  OF  ALTITUDE  ON  THE  UNFOLDING  DIS- 
TANCE AT  CONSTANT  VELOCITY  AND  CONSTAN1 
DYNAMIC  PRESSURE  FOR  n - 1/2  AND  n • 0.632*6. 


FIG  1*  VARIATION  OF  STEADY-STATE  DRAG.  Fr  AND  MAXI- 
MUM OPENING  SHOCK  WITH  ALTITUDE  FOR  CONSTANT 
VELOCITY  AND  CONSTANT  DYNAMIC  PRESSURE 


area  at  the  beginning  of  Inflation.  Once 
t0  has  been  calculated,  a correction  can 
be  applied,  based  upon  what  is  known  about 
the  initial  conditions. 
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Case  A - When  the  Initial  projected 
area  Is  known 


than  the  maximum  extensibility,  e , of 
the  materials.  max 


Development  of  the  analysis  In  the 
elastic  phase  of  inflation  Is  similar  to 
the  technique  used  In  the  unfolding  phase. 
Newton's  second  law  of  motion  Is  used, 
together  with  the  drag  area  ratio  signa- 
ture and  mass  ratio 

CDSo  \V 


corrected 


“calculated 


(15) 


Case  B - When  the  Initial  drag  area 
is  known 

.(1±)6 

CDSo  Vo* 


(16) 

The  mass  ratio  should  row  be  adjusted  for 
the  corrected  t0  before  velocity  and  force 
profiles  are  determined . 

VI.  Opening-Shock  Force,  Velocity  Ratio, 
and  Inflation  Time  During  the 
Elastic  Phase  of  Parachute  Inflation 


The  mass  ratio,  M,  is  an  Important 
parameter  In  parachute  analysis.  For 
values  of  M < < 4/21,  the  maximum  opening- 
shock  force  occurs  early  in  the  inflation 
process,  and  the  elastic  properties  of 
the  canopy  are  not  significant.  As  the 
mass  ratio  approaches  M = 4/21,  the 
magnitude  of  the  opening- shock  force 
Increases,  and  the  time  of  occurrence 
happens  later  In  the  deployment  sequence. 
For  mass  ratios  M > 4/21,  the  maximum 
shock  force  will  occur  after  the  reference 
time,  to.  Parachutes  designed  for  high 
mass  ratio  operation  must  provide  a 
structure  of  sufficient  constructed 
strength,  Fc,  so  that  the  actual  elon- 
gation of  the  canopy  under  load  Is  less 


which  Is  still  valid, as  shown  In  Figure  7 


v 

Integrating  and  solving  for  — 

Vs 


where  — Is  the  velocity  ratio  of  the 
Vs 

unfolding  process  at  time  t = tQ. 


1 + - 


+ ni.1.---!!  + 


(18) 


The  Instantaneous  shock  factor  in  the 
elastic  phase  becomes 


The  end  point  of  the  Inflation  process 
depends  upon  the  applied  loads,  elasticity 
of  the  canopy,  and  the  constructed  strength 
of  the  parachute.  A linear  load  elongation 
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relationship  is  utilized  to  determine  the  The  next  force  in  the  series  at  constant  q 
maximum  drag  area. 

£ _ Cmax  F1  = X0FS  ~ 


r = t1  + eo)£ 


The  force,  F,  is  initially  the  instanta- 
neous force  at  the  end  of  the  unfolding 
process 


F = Xr,Fc, 

o s 


where  XQ  is  the  shock  factor  of  the 
unfolding  phase  at  t = tc 


1+±P 

M L 


+ f 


Since  the  inflated  shape  is  defined,  the 
drag  coefficient  is  considered  to  be 
constant,  and  the  instantaneous  force  is 
proportional  to  the  dynamic  pressure  and 
projected  area.  The  maximum  projected 
area  would  be  developed  if  the  dynamic 
pressure  remained  constant  during  the 
elastic  phase.  Under  very  high  mass 
ratios,  this  is  nearly  the  case  over  this 
very  brief  time  period;  but  as  the  mass 
ratio  decreases,  the  velocity  decay  has 
a more  significant  effect.  The  simplest 
approach  for  all  mass  ratios  is  to 
determine  the  maximum  drag  area  of  the 
canopy  as  if  elastic  Inflation  had 
occurred  at  constant  dynamic  pressure. 

Then  utilizing  the  time  ratio  determined 
as  an  end  point,  intermediate  shock 
factors  can  be  calculated  from  equation 
(19)  and  maximum  force  assessed. 

The  initial  force,  X0FS,  causes  the 
canopy  to  increase  in  projected  area.  The 
new  projected  area  in  turn  increases  the 
total  force  on  the  canopy  which  produces 
a secondary  projected  area  Increase.  The 
resulting  series  of  events  are  resisted 
by  the  parachute  materials.  The  parachute 
must,  therefore,  be  constructed  of  suffi- 
cient strength  to  prevent  the  elongation 
of  the  materials  from  exceeding  the 
maximum  elongation. 


Subsequent  elongations  ir.  the  system  can 
be  shown  to  be 


E1  = co  (1  + Eo> 

C2  - co  + «o  + <o>2>£ 


The  required  canopy  constructed  strength 
can  be  determined  for  a given  set  of 
deployment  conditions.  The  limiting  value 
of  the  series  (e^)  determines  the  end 
point  time  ratio. 


(h\6m5* 
\V  y 


^ = (1  + h)2 


o-fisr-" 


Figure  lb  illustrates  the  maximum  drag 
area  ratio  as  a function  of  e0. 
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FIG.  15  MAXIMUM  DRAG  AREA  RATIO  Vi.  INITIAL 
ELONGATION 
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VII.  Application  of  Cloth  Permeability 
to  the  Calculation  of  the  Inflation 
Time  of  Solid  Cloth  Parachutes 


The  mass  outflow  througn  the  pressur- 
ized region  of  an  inflating  solid  cloth 
parachute  at  any  instant  is  dependent  upon 
the  canopy  area  which  is  subjected  to  air- 
flow and  the  rate  of  airflow  through  that 
area.  The  variation  of  pressurized  area  as 
a function  of  reference  time,  tD,  was 
earlier  assumed  to  be  proportional  to  the 
instantaneous  drag  area  ratio,  leaving  the 
rate- of- airflow  problem  to  solve.  The 
permeability  parameter  of  cloth  was  a 
natural  choice  for  determining  the  rate  of 
airflow  through  the  cloth  as  a function  of 
pressure  differential  across  the  cloth. 
Heretofore,  these  data  have  been  more  of 
a qualitative,  rather  than  quantitative, 
value.  A new  method  of  analysis  was 
developed  wherein  a generalized  curve  of 
the  form  P = k(AP)n  was  fitted  to  cloth 
permeability  data  for  a number  of  different 
cloths  and  gives  surprisingly  good  agree- 
ment over  the  pressure  differential  range 
of  available  data.  The  pressure  differ- 
ential was  then  related  to  the  trajectory 
conditions  to  give  a generalized  expression 
which  can  be  used  in  the  finite  mass  ratio 
range,  as  well  as  the  infinite  mass  case. 
The  permeability  properties  were  trans- 
formed into  a mass  flow  ratio,  M',  which 
shows  agreement  with  the  effective  porosity 
concept. 

Measured  and  calculated  permeability 
pressure  data  for  several  standard  cloths 
are  illustrated  in  Figure  16.  This  method 
has  been  applied  to  various  types  of  cloth 
between  the  extremes  of  a highly  permeable 
3-momme  silk  to  a relatively  impervious 
parachute  pack  container  cloth  with 
reasonably  good  results,  see  Figure  17. 

The  canopy  pressure  coefficient,  C-, 
is  defined  as  the  ratio  of  the  pressure 
differential  across  the  cloth  to  the 
dynamic  pressure  of  the  free  stream. 


C — a Internal)  - P(external) 
p q 1/2  pV2 

(25) 

where  V is  based  on  equation  (7). 

The  permeability  expression,  P = k(aP)n 
becomes 


P 


*(Cp 


(26) 


Although  some  progress  has  been  made 
by  Melzlg  and  others  on  the  measurement  of 
the  variation  of  the  pressure  coefficient 
on  an  actual  Inflating  canopy,  this 
dimension  and  its  variation  with  time  are 
still  dark  areas  at  the  time  of  this 
writing.  At  the  present  time,  a constant 
average  value  of  pressure  coefficient  is 


used  in  these  calculations.  Figure  18 
presents  the  effect  of  pressure  coefficient 
and  altitude  on  the  unfolding  time  for 
constant  deployment  conditions. 

It  is  well  known  that  the  inflation 
time  of  solid  cloth  parachutes  decreases 
as  the  operational  altitude  increases. 

This  effect  can  be  explained  by  consid- 
ering the  ratio  of  the  mass  outflow 
through  a unit  cloth  area  to  the  mass 
inflow  through  a unit  mouth  area. 

,,,  mass  outflow 

M'  - mass  flow  ratio  = 

mass  inflow 


where 

mass  outflow  = P slugs.  (per.  ftg cloth  area) 
" ft  -sec 


and 


mass  inflow  = v (per  ft?  inflow  area) 

P ft  -sec 


REPRODUCED  FROM  REFERENCE  (4) 

FIG.  16  NOMINAL  POROSITY  OF  PARACHUTE  MATERIAL 
VS  DIFFERENTIAL  PRESSURE. 
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Therefore,  the  mass  flow  ratio  becomes 


VO  V 


M' 


(27) 


Effective  porosity,  C,  is  defined  as 
the  ratio  of  the  velocity  through  the 
cloth,  u,  to  a fictitious  theoretical 
velocity,  v,  which  will  produce  the  par- 
ticular iP=  l/2ov2. 


FIG  17  COMPARISON  OF  MEASUREOANO  CALCULATED 

PERMEABILITY  FOR  RELATIVELY  PERMEABLE  AND 
IMPERMEABLE  CLOTHS 


AVERAGE  CANOPY  PRESSURE  COEFFICIENT 
DURING  INFLATION  INCLUDING  THE  VENT 


CP 


FIG.  18  EFFECT  OF  PRESSURE  COEFFICIENT  AND  ALTITUDE 
ON  THE  UNFOLDING  TIME. 


effective  porosity,  C = — (28) 

v 


Comparison  of  the  mass  flow  ratio  and 
previously  published  effective  porosity 
data  is  shorn  in  Figure  19.  The  effects 
of  altitude  and  velocity  on  the  mass  flow 
ratio  are  presented  in  Figures  20,  and  21 
for  constant  velocity  and  constant 
altitude.  The  decrease  of  cloth  perme- 
ability with  altitude  is  evident. 

The  permeability  constants  "k"  and  "i.-' 
can  be  determined  from,  the  permeability 
pressure  differential  data  as  obtained 
from  an  instrument  such  as  a Frazier 
Permeameter.  Two  data  points,  "A"  and 


REPRODUCED  FROM  REFERENCE  (41 


FIG.  19  THE  EFFECTIVE  POROSITY  OF  PARACHUTE 
MATERIALS  VS.  DIFFERENTIAL  PRESSURE 
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I 


"B,  " are  selected  ir.  such  a manner  that 
point  "A"  is  in  a lov.-pressure  zone  below 
the  knee  of  the  curve,  and  point  "B"  is 
located  in  the  upper  end  of  the  high- 
pressure  zone,  as  shown  in  Figure  22. 

The  two  standard  measurements  of  1/2 
inch  of  water  and  20  inches  of  water 
appear  to  be  good  data  points  if  both  are 


FIG.  20  EFFECT  OF  ALTITUDE  ON  MASS  FLOW  RATIO  AT 
CONSTANT  VELOCITY 


FIG  21  EFFECT  OF  VELOCITY  ON  MASS  FLOW  RATIO  AT 
CONSTANT  DENSITY 


FIG.  22  LOCATION  OF  DATA  POINTS  FOR  DETERMINATION 
OF  "k"  AND  "|T. 


available  on  the  same  sample, 
tuting  the  data  from  points  " 

into  P = k(AP)n  : 


n = 


tn 


A" 


Substi- 
and  "B 


II 


(29) 


k = 


(30) 


VIII.  Determination  of  the 
Parachute  Included  Volume 
and  Associated  Air  Mass 


Before  the  reference  time,  tQ,  and 
inflation  time,  tf,  can  be  calculated, 
the  volume  of  atmosphere,  Vo,  which  is 
to  be  collected  during  the  inflation 
process  must  be  accurately  known.  This 
requirement  dictates  that  a realistic 
Inflated  canopy  shape  and  associated 
volume  of  atmosphere  be  determined. 

Figure  23  was  reproduced  from  reference 
(5).  The  technique  of  using  lampblack 
coated  plates  to  determine  the  airflow 
patterns  around  metal  models  of  inflated 
canopy  shapes  was  used  by  the  investigator 
of  reference  (5)  to  study  th»  stability 
characteristics  of  contemporary  parachutes, 
i.e.,  19^3.  A by-product  of  this  study 
is  that  it  is  clearly  shown  that  the 
volume  of  air  within  the  canopy  bulges 
out  of  the  canopy  mouth  (indicated  by 
arrows)  and  extends  aheaa  of  the  canopy 
hem.  This  volume  must  be  collected  during 
the  inflation  process.  Another  neglected, 
but  significant,  source  of  canopy  volume 
exists  ir.  the  billowed  portion  of  the  gore 
panels. 
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- — * 

HEMISPHERE 


VENT  PARACHUTE 

REPRODUCED  FROM  REFERENCE  (5) 


FIG.  73  AIRFLOW  PATTERNS  SHOWING  AIR  VOLUME 
AHEAD  OF  CANOPY  HEM 


The  steady-state  canopy  shape  has  been 
observed  In  wind-tunnel  and  field  tests  to 
be  elliptical  in  profile.  Studies  of  the 
inflated  shape  and  included  volume  of 
several  parachute  types  (flat  circular, 

10  percent  extended  skirt,  elliptical, 
hemispherical,  ring  slot,  ribbon,  and 
cross)  are  documented  in  references  (6) 
and  (7).  These  studies  demonstrated  that 
the  steady-state  profile  shape  of  inflated 
canopies  of  the  various  types  can  be 
approximated  to  be  two  ellipses  of  common 
major  diameter,  2a,  and  dissimilar  minor 
diameters,  b and  b' , as  shown  in  Figure  24. 
It  was  also  shown  that  the  volume  of  the 
ellipsoid  of  revolution  formed  by  revolving 
the  profile  shape  about  the  canopy  axis 
was  a good  approximation  of  the  volume  of 
atmosphere  to  be  collected  during  canopy 
inflation  and  included  the  air  volume 
extended  ahead  of  the  parachute  skirt  hem 
together  with  the  billowed  gore  volume. 


— na 


»t  L CONTOUR 


-h"  \[  EQUIVALENT  BILLOW  VOLUME 

f EXTENDED  AIR  VOLUME 

- b — 

REPRODUCED  FROM  REFERENCE  (6) 

FIG.  24  PARACHUTE  CROSS  SECTION  NOMENCLATURE 

Tables  I and  II  are  summaries  of  test 

results  reproduced  from  references  (6)  and 

(7),  respectively,  for  the  convenience  of 

the  reader. 
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TABLE  I SUMMARY  OF  PARACHUTE  SHAPE  TEST  RESULTS 
FOR  12-GORE  AMD  16-GORE  CONFIGURATIONS 


Parachute 

Type 

No.  of 
Gor  at 

Suspension 
Lma  Length 
inchat 

Valocity 

mph  fpt 

LMA 

B 

Ax  at  Ratio 

b b' 

a S 

b b' 

S S 

Volume  in 

VH  VC 

3 

Vo 

v. 

VH 

Flat  Circular 

12 

50 

73 

645 

650 

PI 

.6115 

8817 

1.4932 

4476 

4481 

6980 

1 56 

16 

50 

73 

663 

669 

5558 

.9039 

1.4597 

4450 

4100 

7325 

1.65 

10%  Extended 

12 

34 

100 

147 

8 63 

852 

881 

6424 

8860 

1.5284 

3928 

4400 

6783 

1.73 

Skirt 

16 

34 

17 

25 

654 

640 

785 

5580 

8502 

1.4082 

4051 

3920 

6197 

1.53 

Elliptical 

12 

34 

75 

110 

916 

812 

5626 

9667 

1.5283 

3322 

5405 

16 

34 

17 

25 

875 

800 

6169 

8163 

1.4332 

2726 

4406 

Hemispherical 

12 

34 

125 

183 

996 

1.0005 

9080 

1 9085 

6224 

8666 

16 

34 

75 

110 

994 

9129 

9380 

1.8509 

5921 

8370 

Rmgsiot 

12 

34 

25 

37 

607 

654 

853 

6566 

8735 

1.530 

3800 

3650 

5903 

1.56 

16%  Gaomatric 

12 

34 

100 

147 

616 

663 

922 

6566 

8735 

1.530 

3800 

4198 

6166 

1 62 

Porosity 

12 

34 

200 

293 

637 

686 

918 

6566 

8735 

1.530 

3800 

4624 

6826 

1.90 

16 

34 

25 

37 

611 

658 

827 

6004 

8890 

1.4894 

3800 

3763 

5685 

1.50 

16 

34 

100 

147 

617 

664 

864 

6004 

8890 

1.4894 

3800 

3985 

6030 

1.59 

16 

34 

200 

293 

645 

695 

844 

6004 

8890 

1 4894 

3800 

4430 

6897 

1.82 

Ribbon 

12 

34 

25 

37 

586 

632 

859 

6558 

8768 

1.5'  ’6 

3800 

3323 

5336 

1.40 

24%  Gaomatric 

12 

34 

100 

147 

615 

663 

837 

6558 

8768 

1.5326 

3800 

3714 

6163 

1 62 

Porosity 

12 

34 

200 

293 

632 

681 

877 

6558 

8768 

1.5326 

3800 

4280 

6683 

1.76 

16 

34 

25 

37 

603 

650 

797 

5570 

8578 

1.4148 

3800 

3438 

5358 

1.41 

18 

34 

100 

147 

626 

674 

791 

5570 

8578 

1.4148 

3800 

3804 

5983 

1.57 

16 

34 

200 

293 

648 

698 

781 

5570 

8578 

1.4148 

3800 

4164 

6656 

1.75 

Croat  Chute 

■ ■ 

34 

25 

37 

710 

543 

1 242 

8867 

1 2776 

2.1643 

1928 

3768 

5798 

3.01 

"/L  ■ 264 

34 

100 

147 

707 

540 

1 270 

8867 

1 277b 

2 1643 

1928 

3810 

5712 

296 

34 

200 

293 

716 

547 

1 285 

8867 

1.2776 

2 1643 

1923 

4212 

5925 

3.07 

47 

25 

37 

759 

580 

1 113 

8494 

1 2512 

2 1006 

1928 

4052 

6868 

356 

47 

100 

147 

729 

557 

1 205 

8494 

1 2512 

2 1006 

1928 

3973 

5958 

309 

47 

200 

293 

775 

592 

1 no 

8494 

1 2512 

2 1006 

1928 

4292 

7303 

3 79 
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TABLE  II  SUMMARY  OF  PARACHUTE  SHAPE  TEST  RESULTS 
FOR  24-GORE  AND  30-GORE  CONFIGURATIONS 


Parachute 

Type 

No  of 

Goret 

Sutpanuon 
Lma  Langth 

inches 

Velocity 
mph  fpt 

Scale  Factor.  K 

2a  2a 

°. 

N 

a 

b 

i 

Axes  Ratio 

k‘  * b 

a s 4 ; 

Volume  in 

V V 

VH  c 

v. 

v. 

VH 

Flat  Circulars 

24 

34 

50 

73 

677 

679 

795 

5758 

8126 

1 3884 

4362 

4695 

7273 

1.67 

30 

34 

17 

25 

668 

669 

827 

6214 

7806 

1 4020 

4342 

4626 

7027 

1 62 

10%  Extended* 

24 

34 

100 

147 

666 

648 

834 

5949 

8771 

1 4720 

4138 

4446 

6930 

1.67 

Skirt 

30 

34 

17 

25 

650 

633 

825 

6255 

7962 

1 4127 

4172 

4076 

6265 

1.60 

Ring  Slot  16% 

24 

34 

25 

37 

663 

665 

824 

5800 

9053 

1 4853 

3591 

3878 

3031 

1 68 

Geometrically 

24 

34 

100 

147 

680 

682 

819 

5800 

9053 

1 4863 

3591 

4079 

6510 

1.81 

Porout 

24 

34 

200 

293 

694 

696 

809 

5800 

9063 

14863 

3591 

4270 

6924 

1.93 

30 

34 

25 

37 

677 

878 

788 

5800 

9053 

1 4863 

3582 

3826 

6404 

179 

30 

34 

100 

147 

684 

685 

802 

5800 

9053 

1 4853 

3582 

4023 

6588 

1 84 

30 

34 

200 

293 

698 

694 

800 

5800 

9053 

1 4863 

3582 

4260 

7012 

1 96 

Ribbon  24% 

24 

34 

25 

37 

671 

673 

770 

5960 

8187 

1.4167 

3591 

3591 

5968 

1 66 

Geometrically 

24 

34 

100 

147 

676 

678 

813 

5960 

8187 

1.4167 

3591 

3927 

6097 

1.70 

Porout 

24 

34 

200 

293 

687 

689 

804 

5960 

8187 

1.4167 

3591 

4061 

6389 

1 78 

30 

25 

37 

655 

657 

782 

6021 

8463 

1 4484 

3682 

3396 

5666 

1 58 

30 

100 

147 

669 

670 

764 

6021 

8463 

1.4464 

3562 

3622 

6022 

1 68 

30 

H3I 

200 

293 

677 

679 

823 

6021 

8463 

1 4484 

3682 

4002 

6256 

1 75 

*tmce  this  parachute  was  'breathing''  during  the  test.  several  photographs  wart  taken  at  each  speed  The  data  were  reduced  from  the  photograph  which  mot t 
reasonably  appeared  to  rapraaant  the  epuilibruim  itata 
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A 


c 


\t 

*Vo 


a 

73 


b 


b' 


F 

F„ 


K 

k 


m 

M 


M' 
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List  of  Symbols 

- Steady-state  pro jected^area  of  the 
Inflated  parachute,  ft- 

- Instantaneous  canopy  mouth  area,  ft^ 

- Steady-state  Inflated  mouth  atoa, 
ft? 

- Acceleration,  ft/sec^ 

- Maximum  inflated  parachute  diameter 
of  gore  mains earn,  ft 

- Minor  axis  of  the  ellipse  bounded 
by  the  major  axis  (2a)  and  the  vent 
of  the  canopy,  ft 

- Minor  axis  of  the  ellipse  which 
includes  the  skirt  hem  of  the 
canopy,  ft 

- Effective  porosity 

- Parachute  coefficient  of  drag 

- Parachute  pressure  coefficient, 
relates  internal  and  external 
pressure  (aP)  on  canopy  surface  to 
the  dynamic  pressure  of  the  free 
stream 

- Nominal  diameter  of  the  aerodynamic 

decelerator  = , ft 

- Instantaneous  force,  lbs 

- Steady-state  drag  force  that  would 
be  produced  by  a fully  cpen  para- 
chute at  velocity  Vs,  lbs 

- Constructed  strength  of  the  para- 
chute, lbs 

- Maximum  opening-shock  force,  lbs 

O 

- Gravitational  acceleration,  ft/sec 

- Permeability  constant  of  canopy 
cloth 

- Mass,  slugs 

- Mass  ratio  - ratio  of  the  mass  of 
the  retarded  hardware  (including 
parachute)  to  a mass  of  atmosphere 
contained  In  a right  circular 
cylinder  of  length  (VstQ),  face  area 

(cdso)»  and  density  (°) 

- Mass  flow  ratio  - ratio  of  atmo- 
sphere flowing  through  a unit  cloth 
area  to  the  atmosphere  flowing 
tirough  a unit  inlet  area  at 
arbitrary  pressure 

- Permeability  constant  of  canopy 
cloth 


78-189 

P - Cloth  permeability  - rate  ol  air- 
flow through  a cloth  at  an  arbi- 
trary differential  pressure, 
ft 3/ft ? sec 

p 

q - Dynamic  pressure.  It/ ft' 

S - Instantaneous  inflated  canopy 

surface  area,  ft2 

p 

S0  ' - Canopy  surface  area,  ft 

t - Instantaneous  time,  sec 

t - Reference  time  when  the  parachute 

has  reached  the  design  drag  area 
for  the  first  time,  sec 

t , - Canopy  inflation  time  when  the 

Inflated  canopy  has  reached  its 
maximum  physical  size,  sec 

u - Air  velocity  through  cloth  in 
effective  porosity,  ft/sec 

v - Fictitious  theoretical  velocity 

used  in  effective  porosity,  ft/sec 

V - Instantaneous  system  velocity,  ft/sec 

V„  - System  velocity  at  the  time  t = t , 

° ft/sec  0 

V - System  velocity  at  the  end  of 

suspension  line  stretch,  ft/sec 

V - "olume  of  air  which  must  be  collec- 

ted during  the  inflation  process, 

ft3 

W - Hardware  weight,  lb 

x^  - Instantaneous  shock  factor 

X - Shock  factor  at  the  time  t = t 

o o 

p - Air  density,  slugs/ft3 

p - Ratio  of  parachute  projected  mouth 

area  at  line  stretch  to  the  steady- 
state  projected  area 

e - Instantaneous  elongation 

e - Maximum  elongation 

max 

£ - Initial  elongation  at  the  beginning 

of  the  elastic  phase  of  inflation 

S.F.  - Parachute  safety  factor  = F /F^ 
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Appendix  C 

A GUIDE  FOR  THE  USE  OF  APPENDIX  B 


At  first  reading,  Appendix  B may  appear  to  be  a complicated  sys- 
tem of  analysis  because  of  the  many  formulas  presented.  Actually, 
once  understood,  the  technique  is  straightforward  and  uncomplicated. 
The  author  has  attempted  to  simplify  the  algebra  wherever  possible. 
This  appendix  presents,  in  semi-outline  form,  a guide  to  the 
sequence  of  calculations  because  the  analysis  does  require  use  of 
formulas  from  the  text,  not  necessarily  in  the  order  in  which  they 
were  presented.  Also,  the  user  can  be  referenced  to  graphs  of  per- 
formance to  illustrate  effects. 

In  order  to  compute  tQ,  other  parameters  must  be  obtained  from 
various  sources. 

I.  Determine  System  Parameters 

1.  C[)Sq,  drag  area,,  ft^  obtained  from  design  requirement. 

2.  Vs,  fps,  velocity  of  system  at  suspension  line  stretch. 

3.  p,  slugs/ft-^,  air  density  at  deployment  altitude. 

4.  W , lb,  system  weight  (including  weight  of  the  parachute) 
from  design  requirements. 

5.  VQ,  ft^,  this  volume  of  air,  which  is  to  be  collected 
during  inflation,  is  calculated  from  the  steady-state  inflated 
shape  geometry  of  the  particular  parachute  type.  The  nomenclature 
is  described  in  Figure  24,p.B-13.  When  D0  or  Dp  is  known,  a can  be 
calculated  from  data  in  Table  I and  Table  II,  p.  B-14,  lor  various 
parachute  types  and  number  of  gores.  Then  the  geometric  volume  VQ 
can  be  calculated  by  Equation  (31),  p.  B-13,  with  appropriate 
values  of  b/a  and  b'/a  from  the  tables. 

6.  Amo»  ft 2,  steady-state  canopy  mouth  area 
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where  n/a,  b/a,  and  b'/a  are  available  from  Tables  I and  II  for  the 
particular  type  of  parachute  and  number  of  gores. 

7.  Aso,  ft2,  canopy  surface  area  = 

8.  C , pressure  coefficient,  see  Figure  18,  p.  B-ll.  A 
constant  Cp  = 1.7  for  all  altitudes  seems  to  yield  acceptaole 
results . 

9.  Constants  k and  n are  derived  from  measurements  of  the  air 
flow  througn  the  cloth.  Only  k is  needed  for  Equation  (14),  but  n 
is  also  required  for  Equation  (13) . These  parameters  can  be  deter- 
mined for  any  cloth  using  the  technique  described  beginning  on 

p.  B-ll.  The  two-point  method  is  adequate  if  the  AP  across  the 
cloth  is  in  the  range  of  AP  for  actual  operation.  Check-points  of 
cloth  permeability  can  be  measured  and  compared  to  calculated 
values  to  verify  agreement.  If  the  data  is  to  be  extrapolated  to 
operational  AP ' s greater  than  measured,  a better  method  of  deter- 
mining k and  n from  the  test  data  would  be  a least  squares  fit 
through  many  data  points.  This  way  errors  due  to  reading  either 
of  the  two  points  are  minimized. 

11.  Step  1 

Calculate  the  reference  time  tQ  by  use  of  Equations  (13)  or 
(14),  p.  B-6.  If  the  deployment  altitude  is  50,000  feet  or  higher, 
Equation  (14)  is  preferred  due  to  its  simplicity.  For  altitudes 
from  sea  level  to  50,000  feet.  Equation  (13)  is  preferred.  Figure 

12,  p.  B-7,  shows  the  effect  of  altitude  on  tQ  and  can  be  taken  as 

a guide  for  the  user  to  decide  whether  to  use  Equation  (13)  or  (14). 
One  should  keep  in  mind  that  the  opening  shock  force  can  be  a 
strong  function  of  inflation  time,  so  be  as  realistic  as  possible. 

If  Equation  (13)  is  elected,  the  method  in  use  at  the  NSWC/WO  is 
to  program  Equation  (13)  to  compute  the  parachute  volume,  V^,,  for 
an  assumed  value  of  tQ.  Equation  (14),  because  of  its  simplicity, 
can  be  used  for  a first  estimate  of  tQ  at  all  altitudes.  The  com- 
puted canopy  volume  is  then  compared  to  the  canopy  volume  calculated 
from  the  geometry  of  the  parachute  as  per  Equation  (31),  p.  B-  3. 

If  the  volume  computed  from  the  mass  flow  is  within  the  volume 
computed  from  the  geometry  within  plus  or  minus  a specified  delta 
volume,  the  time  t0  is  printed  out.  If  not  within  the  specified 
limits,  tQ  is  adjusted,  and  a new  volume  calculated.  For  a 35-foot 
D0 , T-10  type  canopy,  I use  plus  or  minus  10  cubic  feet  in  the 

volume  comparison.  The  limit  would  be  reduced  for  a parachute  of 

smaller  D&. 

If  V0  calculated  = VD  geometry  ± 10,  then  print  answer. 

If  VQ  calculated  f VQ  geometry  + 10,  then  correct  tQ  as  follows 

~Y0  geometry 

X°  ' Y„  calculated  (C-2) 
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The  new  value  of  tQ  is  substituted  in  the  "do  loop"  and  the 
volume  recomputed.  This  calculation  continues  until  the  required 
volume  is  within  the  specified  limits. 

III.  Calculate  tD  corrected  for  initial  area.  The  tp  of  Section 
II  assumes  that  the  parachute  inflated  from  a zero  initial  are?. 

If  this  is  a reasonable  assumption  for  the  particular  system  under 
study,  then  the  mass  ratio  can  be  determined  from  Equation  (6) , 

p.  B-3.  For  n = 0 if  the  value  of  M l 0.19,  then  a finite  state 
of  deployment  exists,  and  the  time  ratio  of  occurrence  and  the 
maximum  shock  factor  can  be  determined  from  Equations  (9)  and  (10) 
respectively  on  p.  B-4.  If  n ^ 0,  then  the  limiting  mass  ratio 
for  finite  operations  will  rise  slightly  as  described  in  Appendix 
D.  Figures  D-l  and  D-2  illustrate  the  effects  of  initial  area  on 
limiting  mass  ratios  and  shock  factors  respectively.  If  the  mass 
ratio  is  greater  th  a n the  limiting  mass  ratio  (ML) , then  the 
maximum  shock  force  occurs  at  a time  greater  than  t0  and  the 
elasticity  of  the  materials  must  be  considered  (see  Section  IV) . 

If  n ^ 0,  then  the  reference  time  tQ  will  be  reduced,  and 
the  mass  ratio  will  rise  due  to  partial  inflation  at  the  line 
stretch.  Figures  9 and  10,  p.  B-5,  illustrate  the  effects  of 
initial  area  on  the  velocities  and  shock  factor  during  the 
"unfolding"  inflation.  Equation  (15),  p.  B-8,  can  be  used  to 
correct  tQ  calculated  for  where  n = Ai/Ac.  If  the  initial  value 
of  drag  area  is  known,  Equation  (16),  p.  B-8,  can  be  used  to 
correct  t and  rechecked  for  limiting  mass  ratios  versus  n in 
Appendix  D. 

IV.  Opening  shock  calculations  in  the  elastic  phase  of  inflation. 
It  has  been  considered  that  from  time  t = 0 to  t = tQ  the  para- 
chute has  been  inelastic.  At  the  time  t = tQ  the  applied  aero- 
dynamic load  causes  the  materials  to  stretch  and  the  parachute 
canopy  increases  in  size.  The  increased  size  results  in  an 
increase  in  load,  which  causes  further  growth,  etc.  This  sequence 
of  events  continues  until  the  applied  forces  have  been  balanced  by 
the  strength  of  materials.  The  designer  must  insure  that  the 
constructed  strength  of  the  materials  is  sufficient  to  resist  the 
applied  loads  for  the  material  elongation  expected.  Use  of 
materials  of  low  elongation  should  result  in  lower  opening  shock 
forces  as  CDSmax  is  reduced. 

When  the  mass  ratio  of  the  system  is  greater  than  the  limiting 
mass  ratio,  the  elasticity  of  the  materials  and  material  strength 
determine  the  maximum  opening  shock  force.  The  maximum  elongation 
emax  anc^  the  ultimate  strength  of  the  materials  are  known  from 
tests  or  specifications.  The  technique  begins  on  p.  B-8. 

At  the  time  t = tQ  calculate  the  following  quantities  for  the 
particular  values  of  M and  n- 

a.  V0/Vs  from  Equation  (18),  p.  B-8. 
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b.  XQ  from  Equation  (22),  p.  B-9. 

c.  eQ  from  Equation  (23),  p.  B-9. 

d.  Determine  CDSmax/CDS0  from  Figure  15,  p.  B-9. 

e.  Calculate  the  inflation  time  ratio  tf/tQ  from  Equation 
(24),  p.  B-9. 

f.  Calculate  the  maximum  shock  factor  from  Equation  (19), 
p.  B-8. 

g.  Calculate  the  opening  shock  force  “ \FS  where 

FS  - 

h.  Calculate  filling  time,  tf(sec) 

V.  In  order  to  simplify  the  required  effort,  the  work  sheets  of  Table 
C — 1 are  included  on  pages  C-5  through  C-9  to  aid  the  engineer  in 
systematizing  the  analysis.  The  work  sheets  should  be  reproduced  to 
provide  additional  copies. 
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Table  C-l.  Opening  Shock  Force 

CALCULATION  WORK  SHEETS  \ JV 


1.  Parachute  type  - 

2.  System  parameters 

a.  System  weight,  W (lb) 

b.  Gravity  , g (ft/sec2) 

c.  Deployment  altitude  (ft) 

d.  Deployment  air  density,  p (slugs/ft3) 

e.  Velocity  at  line  stretch,  Vs  (fps) 

f.  Steady  state  canopy  data 


2a 


(1) 

Diameter,  DQ  (ft) 

(2) 

Inflated  diameter,  2a 

(ft); 

(3) 

Surface  area,  SQ  (ft2) 

• - D 

’ 4 

(4) 

Drag  area,  CDS0  (f t2) ; 

CD 

(5) 

Mouth  area,  A^q  (ft) 

amo  ■ nt 


(6)  Volume,  VQ*  (ft3) 

[t'f] 

g.  Cloth  data 

(1)  K J Calculate  using  technique  beginning  on 

(2)  n ) p.  B-ll.  Note:  Permeability  is  usually 

measured  as  ft3/ft2/min. 

For  these  calculations 
permeability  must  be 
expressed  as  ft3/ft2/sec 


* Data  for  these  calculations  are  listed  in  Tables  1 
and  2,  p.  B-14. 
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W 

s 

p 

V. 


2i 

So 

sDs0 


amo 


lb. 

ft/wc2 

ft. 

•lujt/ft2 

fp». 

ft. 

ft. 

ft. 

ft.2 


ft.2 


ft.J 
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Table  C-l.  Opening  Shock  Force 

(cont'd) 

(3)  determine  maximum  elongations  from  ^ 

pull  test  data  of  joints,  seams,  lines,  etc.  Use 

minimum  ^ determined  from  tests. 

(4)  Cp;  pressure  coefficient 

h.  Steady  state  drag,  Fs  (lb)  “ Fs  - ~pV|CQSo 


i.  Parachute  constructed  strength,  Fc  (lb);  deter- 
mined from  data  on  efficiency  of  seams,  joints,  lines. 
Constructed  strength  is  the  minimum  load  required  to 
fail  a member  times  the  number  of  members. 

3.  Force  calculations 


a.  Calculate  tQ  for  n = 0;  eq.  14,  p.  B-6. 


|-  paYo  r D 0 

I 2W  L AMo  ' *50  * ^CpP 


S’)*  ] 

2 ' -1 


Check  Figure  12,  p.  B-7  for  advisability  of 
using  eq.  13,  p.  B-6. 

b.  If  n = 0,  proceed  with  steps  c through  e. 

If  n i 0,  go  to  step  f. 

c.  Mass  ratio,  M;  eq.  6,  p.  B-3 


*«V.CDS. 


d.  If  M < 4/21  for  n » 0,  then  finite  mass 
deployment  is  indicated. 

(1)  Time  ratio  at  Xj  eq.  9,  p.  B-4 


(2)  Max  shock  factor,  xt;  eq. , 10,  p.  B-4 
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Table  C-l.  Opening  Shock  Force 

(Cont'd) 

(3)  Max  shock  force,  Fmax  (lb) 

F - X,  F. 

1 max.  * 

e.  If  M > 4/21;  then  intermediate  mass  or  infinite 
mass  deployment  is  indicated  and  the  elasticity  of 
materials  is  involved.  Calculate  the  trajectory  con- 
ditions at  time  t - tQ. 

(1)  Velocity  ratio  @ t - tQ  for  n « 0 


Vs  . . 1 

8 1 + m 


(2)  Shock  factor  XQ  @ t ■ tQ  for  n “ 0 


[1**]  V V 

(3)  Initial  elongation,  eQ;  eq.  23,  p.  B-9 


XoFS 

fo  ’ F~  tmm 
e 


p.  B-9 


(4) 

Determine 

D max  from  pigure 

15,  p.  B-9 

C„S 

Do 

(5) 

Calculate 

inflation  time  ratio. 

tf  ; eq.  24 

t, 


(6)  Calculate  maximum  shock  factor,  x^  eq. 

19,  p.  B-8  /t,\* 

\w 


x,m»«  r V, 


(7)  Calculate  maximum  shock  force,  (lb), 


I 


NSWC/WOL  TR  78-189 


Table  C-l.  Opening  Shock  Force 

„ (cont’d) 


(8)  Inflation  time,  sec  =■  *f  “ *o 


ft) 


f.  If  n I1  0,  correct  tQ  for  initial  area  effects; 

eq.  16,  p.  B-8  [ /CDSj\1'«1 

0 ' | (cs  ) 10 

g.  Mass  Ratio,  M,  eq.  6,  p.  B-3  ^ \ D o / J 


%% 


PJVoVo 


h.  Calculate  limiting  mass  ratio,  M^;  eq.  D-2, 
p.  D-2 


1 f 9 3 3 1 

Mt  ‘ 3(1  -q)‘  14  n * 14  ^ 7 


If  M < , finite  mass  deployment  is  indicated  and  x^  max 

can  be  determined  by  eq.  8,  p.  B-4  by  assuming  values 
of  t/tQ  and  plotting  the  data  using  the  methods  of 
Appendix  D. 

i.  If  M > M^,  then  intermediate  mass  or 
infinite  mass  deployment  is  indicated  and  the 
elasticity  of  materials  is  involved.  Calculate 
the  trajectory  conditions  at  time  t = t0. 


p.  B-8 


(1)  Velocity  ratio  @ t = tQ  for  n t 0 ; eq.  18, 


22,  p.  B-9 


w — _ 

ri  + irHzJ2l!  + a<Ljz» +r)J 

L M L 7 2 


(2)  Shock  factor  XQ  @ t » tQ  for  r|  f 0;  eq. 


1 fcV 

j^.I  J W 


(3)  Initial  elongation,  eQ;  eq.  23,  p.  B-9 


V — »•« 


c s 

(4)  Determine  — from  Figure  15,  p.  B-9 


CDsm«x 

cD*o 
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Table  C-l.  Opening  Shock  Force  (Con 

(5)  Calculate  Inflation  time  ratio,  — ; eq . 24, 

P-  B-9  . i to 


*,  /CDSm„\‘ 


(6)  Calculate  maximum  shock  factor,  x^  max; 
eq.  19,  p.  B-8  / t \ • 


mm  T V, 


(*r 


(7)  Calculate  maximum  shock  force,  Fmflx(lb) 


F~>  ‘ r‘ 


(8)  Calculate  inflation  time,  tf(sec) 
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Appendix  D 

EFFECT  OF  INITIAL  AREA  RATIO  ON  THE  LIMITING  MASS  RATIO 
AND  SHOCK  FACTOR  FOR  THE  FINITE  STATE  OF 
SOLID  CLOTH  PARACHUTE  DEPLOYMENT 

Very  low  mass  ratios  are  indicative  of  finite  mass  parachute 
deployment,  wherein  the  maximum  shock  force  occurs  before  the 
parachute  is  fully  inflated  during  the  unfolding  phase  of  deploy- 
ment. As  the  mass  ratio  is  increased,  the  maximum  shock  force 
occurs  later  in  the  inflation  process.  At  some  value  of  mass 
ratio,  the  maximum  shock  force  will  occur  at  the  time  tQ.  This 
particular  mass  ratio  is  defined  as  the  limiting  mass  ratio  (ML) 
for  finite  mass  deployment.  A further  increase  in  mass  ratio 
will  result  in  the  maximum  shock  force  occurring  after  the  para- 
chute has  achieved  the  design  drag  area  (CdSq)  for  the  first  time. 

Equation  (8),  from  p.  4,  Appendix  B,  defines  the  instantaneous 
shock  factor  during  the  unfolding  phase  of  parachute  deployment. 


This  expression  is  to  be  analyzed  for  the  following  purposes: 

a.  Determine  the  effect  of  the  initial  area  ratio  (n)  on  the 
limiting  mass  ratio  (ML) . 

b.  Determine  the  variations  of  the  instantaneous  shock  factor 
during  the  unfolding  phase  of  deployment  as  a function  of  limiting 
mass  ratio  and  n. 

c.  Determine  the  expression  for  the  time  of  occurrence  of  the 
maximum  shock  force  for  finite  mass  ratios  less  than  the  limiting 
mass  ratio. 

The  maximum  shock  force  occurs  at  the  point  in  finite  mass 
deployment  where  dx^/dt  = 0.  Setting  the  derivative  dxj/dt  = 0 
and  solving  for  mass  ratio  as  a function  of  n and  t/tQ  results  in 
the  following  equality: 


D-1 


I 
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<1- 

6 

| + 2nd  -i?>  | 

8H 

3 

* + ijd-' 

-6H 

>?(1  -r,) 
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ion 


Since  the  limiting  mass  ratio  occurs  at  t/tQ  = 1,  Equation  (D-l) 
can  be  reduced  to: 


M 


L 


1 

3 (1  - T7» 


2. 

14 


V 


(D-2) 


The  effects  of  initial  area  ratio  on  the  limiting  mass  ratio  are 
described  in  Equation  (D-2)  and  Figure  D-l.  Note  that  the  time  of 
occurrence  of  the  maximum  shock  force  for  n = 0 in  Equation  (D-l) 
is : 


max 


(D-3) 


which  is  the  same  as  Equation  (9),  p.  4,  Appendix  B. 

The  variation  of  the  instantaneous  shock  factor  during  the 
unfolding  phase  of  deployment  for  limiting  mass  ratios  is  presented 
in  Figure  D-2.  Initial  area  at  the  beginning  of  inflation  causes 
the  initial  force  to  increase,  but  this  is  compensated  for  by 
reduced  maximum  shock  forces.  As  n increases,  the  initial  loads 
can  be  greater  than  the  maximum  shock  force.  However,  values  of 
n are  usually  snail  and  depend  on  the  deployment  systems  for 
magnitude  and  repeatability.  Values  of  n = 0.4  are  more  repre- 
sentative of  a parachute  being  disreefed  rather  than  initially 
deployed.  This  does  demonstrate,  however,  that  the  analysis  pre- 
sented in  Appendix  B can  be  adapted  to  the  disreefing  of  solid  cloth 
parachutes  by  considering  the  next  stage  to  be  a deployment  with  a 
large  value  of  n.  Variation  in  initial  area  is  one  of  the  causes 
of  variation  in  opening  shock  forces.  The  variation  of  opening 
shock  forces  for  finite  mass,  intermediate  mass,  and  infinite  mass 
states  of  deployment  can  be  evaluated  by  successive  calculations 
with  various  expected  values  of  n. 

For  known  mass  ratios  less  than  the  limiting  mass  ratio,  the 
time  of  occurrence  of  the  maximum  shock  force  can  be  ascertained 
from  Equation  D-l.  If  n approaches  zero,  then  the  time  ratio  of 


D-2 


0 0.1  0.2  0.3  0.4 


INITIAL  AREA  RATIO,  tj 

FIGURE  D-1  EFFECT  OF  INITIAL  AREA  RATIO  ON  THE  LIMITING 
MASS  RATIO  FOR  TKH  FINITE  STATE  OF  PARACHUTE 
DEPLOYMENT  FOR  SOLID  CLOTH  PARACHUTES 


D-3 


INSTANTANEOUS  SHOCK  FACTOR 
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FIGURE  D-2  VARIATION  OF  THE  FINITE  MASS  SHOCK  FACTOR 
DURING  THE  UNFOLDING  PHASE  OF  SOLID  CLOTH 
PARACHUTES  FOR  LIMITING  MASS  RATIOS  AND 
INITIAL  AREA  EFFECTS 


D-4 
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occurrence  of  the  maximum  shock  force  can  be  initially  estimated 
from  Equation  (D-3) , or  determined  by  plotting 

Xi  = f h,  n,  i — 

zo 

as  in  Figure  D-2. 


D-5 
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SYMBOLS 


Steady-state  projected  area  of  the  inflated  parachute,  ft2 


Projected  area  at  line  stretch,  ft2 


Instantaneous  canopy  mouth  area,  ft2 
Steady-state  inflated  mouth  area,  ft^ 


Acceleration,  ft/sec2 


Maximum  inflated  parachute  diameter  of  gore  mainseam,  ft 


Minor  axis  of  the  ellipse  bounded  by  the  major  axis  (2a)  and 
the  vent  of  the  canopy,  ft 


Minor  axis  of  the  ellipse  which  includes  the  skirt  hem  of 
the  canopy,  ft 


Speed  of  sound  at  analysis  altitude,  fps 


CDS0 


Effective  porosity 


Parachute  coefficient  of  drag 


Steady  state  parachute  drag  area,  ft2 


Parachute  pressure  coefficient,  relates  internal  and  external 
pressure  (AP)  on  canopy  surface  to  the  dynamic  pressure  of 
the  free  stream 


S-l 
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Da  - Aerodynamic  diameter,  ft 

Dq  - Nominal  diameter  of  the  aerodynamic  decelerator  = 
/4Sq/tt,  ft 


Fr  - Froude  number 

F - Instantaneous  force,  lbs 


F 

rmax 

g 

k 

m 


Steady-state  drag  force  that  would  be  produced  by  a fully 
open  parachute  at  velocity  Vg,  lbs 

Constructed  strength  of  the  parachute,  lbs 

Maximum  opening-shock  force,  lbs 

Gravitational  acceleration,  ft/sec 

Permeability  constant  of  canopy  cloth 

Mass,  slugs 


M - Mass  ratio  - ratio  of  the  mass  of  the  retarded  hardware 

(including  parachute)  to  a mass  of  atmosphere  contained  in 
a right  circular  cylinder  of  length  (VstQ) , face  area 

(CDS0) / and  density  (p) 

M'  - Mass  flow  ratio  - ratio  of  atmosphere  flowing  through  a 
unit  cloth  area  to  the  atmosphere  flowing  through  a unit 
inlet  area  at  arbitrary  pressure 

Mn  - Mach  number 

n - Permeability  constant  of  canopy  cloth 


P - Cloth  permeability  - rate  of  airflow  through  a cloth  at  an 

arbitrary  differential  pressure,  ft^/ft^  sec 


S-2 
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Dynamic  pressure,  lb/ft2 
Reynolds  number 

Instantaneous  inflated  canopy  surface  area,  ft 
Canopy  surface  area,  ft2 

Instantaneous  time,  sec 

Reference  time  when  the  parachute  has  reached  the  design 
drag  area  for  the  first  time,  sec 

Canopy  inflation  time  when  the  inflated  canopy  has  reached 
its  maximum  physical  size,  sec 

Air  velocity  through  cloth  in  effective  porosity,  ft/sec 

Fictitious  theoretical  velocity  used  in  effective  porosity, 
ft/sec 

Instantaneous  system  velocity,  ft/sec 

System  velocity  at  the  time  t = tQ,  ft/sec 

System  velocity  at  the  end  of  suspension  line  stretch, 
ft/sec 

Volume  of  air  which  must  be  collected  during  the  inflation 
process,  ft^ 

Hardware  weight,  lb 

Instantaneous  shock  factor 

Shock  factor  at  the  time  t = tG 

Air  density,  slugs/ft^ 


S-3 
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Ratio  of  parachute  projected  mouth  are*  at  line  stretch  to 
the  steady-state  projected  area 


Instantaneous  elongation 


-max 


Maximum  elongation 


e0  - Initial  elongation  at  the  beginning  of  the  elastic  phase  of 
inflation 


S.F.  - Parachute  safety  factor  = Fc/F 


c'  max 


„ . . lb-sec 

U - Dynamic  viscosity,  — - — ^ — 
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GLOSSARY  OF  TERMS 


Aerodynamic  diameter.  Preference  length  used  in  the  computation  of 
Reynold's  number  and  Froude  number. 

°A  = 


Aerodynamic  size.  Name  for  steady-state  parachute  drag  area,  CDS0  ■ 


Ballistic  coefficient.  Ratio  of  retarded  weight  to  parachute  steady- 
state  drag  area,  W/CdSq. 


Elastic  phase  of  inflation.  Inflation  time  from  t0<t£tf  when  the 
elasticity  of  parachute  materials  are  resisting  opening  shock  force. 


Finite  mass  state  of  operation.  Maximum  shock  force  occurs  before 
complete  inflation  at  time  t=tQ. 


Intermediate  mass  state  of  operation.  Maximum  shock  force  occurs 
after  the  parachute  has  opened  to  its  steady-state  drag  area  for  the 
first  time.  The  mass  ratio  M^<M<10.  At  low  mass  ratios  maximum 
force  can  be  much  lower  than  infinite  mass  calculations. 


Limiting  mass  ratio.  The  highest  mass  ratio  for  finite  mass  state  of 
operation.  Maximum  shock  occurs  at  t=tQ. 

State  of  parachute  operation.  Relates  to  time  of  occurrence  of  the 
opening  shock  force  and  is  a function  of  / \ . 

Unfolding  phase  of  inflation.  Inflation  time  from  0<t^tQ  when  the 
canopy  is  inflating. 
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